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J> . Abstract 

I An asymptotic theory is developed for general non-integrable boundary quantum field 

theory in 1+1 dimensions based on the Langrangean description. Reflection matrices are 
defined to connect asymptotic states and are shown to be related to the Green functions 
via the boundary reduction formula derived. The definition of the -R-matrix for integrable 
theories due to Ghoshal and Zamolodchikov and the one used in the perturbative approaches 
are shown to be related. 



Introduction 

Two dimensional boundary quantum field theories have been analysed from two different points 
^ . of view, the bootstrap and the perturbative, respectively. 

The former was initiated by Ghoshal and Zamolodchikov in (l| and can be applied to inte- 
grable theories. In such theories there is an infinite number of conserved quantities, which give 
severe restrictions on the allowed physical processes: Besides the usual constraints such as factor- 
ization and purely elastic bulk scattering there is also factorization and purely elastic reflection 
on the boundary. The scattering theory developed in (l| is analogous to axiomatic scattering 
theory (15|: in the in state the particles travel towards the boundary with decreasingly ordered 
momenta, while in the out state, where all the scatterings and reflections have been terminated, 
they travel away from the boundary with decreasingly ordered momenta again. The i?-matrix 
which connects the in and out states is the composition of the individual reflection and the pair- 
wise scattering matrices. The one particle reflection matrices have to obey unitarity, boundary 
Yang-Baxter and boundary crossing relations. Using these relations together with the bootstrap 
condition (0,@]) the model can be solved modulo CDD type ambiguities. We emphasise that in 
the bulk case the axioms of the scattering theory such as unitarity, crossing symmetry |l5| were 
motivated by relativistic field theoretic results based on the perturbative, Lagrangean descrip- 
tion. In the boundary case, however, to our knowledge, no such background is available. In 
the non-linear Schrodinger model with linear boundary condition on the half line is considered 
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and the assumptions of the axiomatic scattering theory are rigorously checked. This model is, 
however, non-relativistic and integrable. 

The perturbative approach to boundary quantum field theories can be applied without the 
assumption of integrability. It was started with the analysis of bulk perturbation [f|, |5| with 
Neumann boundary condition, then extended to boundary perturbations in ||-||12|]. Most of 
these papers deal with comparing exact results, obtained in the aforementioned way for the 
reflection matrices in integrable theories, on one hand, and perturbative results on the other. 
They defined the reflection matrix R(k) through the asymptotic behaviour of the two point 
function of the field $ - creating the particles - far away from the boundary: 

/J, . p -iw(t-t') 
——-(e^MI*-* 7 ! + R(k)e- ik ^ )ix+x,) ) + ... (1) 
£71 2k{UJ) 

No detailed justification has been given for this quantity being the same as the axiomatic R- 
matrix. 

In this paper we develop the asymptotic theory for a scalar field with the most general bulk 
and boundary self-interaction. We derive the boundary reduction formula, which connects the 
reflection matrix to the two point function. As a consequence we can the fill the gap mentioned 
above, that is we are able to connect the .R-matrix of the axiomatic theory to (|l|). 

The paper is organized as follows: We apply the canonical quantization procedure to the free 
theory, in which case the boundary condition is Neumann. The interacting theory is defined 
by means of the adiabatic hypothesis. Asymptotic states and reflection matrices are introduced 
and the simplest physical process of one incoming particle is demonstrated. As the main result 
we derive the boundary reduction formula. Having developed the boundary perturbation theory 
we are able to connect the earlier definitions of the .R-matrix, and finally we conclude on their 
equivalence. The brief explanation of the Feynman rules is placed in Appendix A, while the 
structure of the two point function is analysed in Appendix B. 



The free theory 

The system we are dealing with contains a real scalar field t), living on the half space x < 0. 
The bulk and boundary interactions are described by the action: 



S 



dt 



dx 



-((d t $) 2 -(^$) 2 -m 2 $ 2 )-V($) 



dtU($(0,t)) 



(2) 



The free theory can be obtained by switching off the bulk and the boundary interactions: V(<&) = 
U(Q>) = 0. The equation of motion is the usual bulk free equation, the boundary condition is, 
however, the Neumann one: 

(□ + m 2 )$(x,t) =0 ; d x ${x,t)\ x=o = . 

In solving these equations by Fourier transformation we have to use the complete system of 
functions with this boundary condition: 



dk 
2^ 



cos(fcx)$(M) ; $(M) = $(-M) • 
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The conjugate momentum also satisfies Neumann boundary condition so the canonical commu- 
tation relation reads as follows: 

t), n(x', t)} = iSjsr(x, x) = i5(x — x) + iS(x + x) . 

The creation and annihilation operators which diagonalise the Hamiltonian are 

a(k, t) — ifl(k, t) + u(k)®(k, t) ; a + (k, t) = -if[(k, t) + u(k)$(k, t) , 

where u(k) = \Jk 2 + m 2 . Their commutation relations are 

[a{k, t), a + {k', t)} = 27i2u{k) {5{k - k') + 5{k + k')) . 

Normal ordering is defined as usual: creation operators a + (k,t) are to the left of annihilation 
operators a(k',t). Since the normal ordered Hamiltonian is 

1 f°° ~ dk 
H=- dkuj(k)a + (kA)a(kA) ; dk = — - 

2 loo K J K J K } 2n2u(k) 

the time dependence can be determined exactly: a + (k, t) = a + (fc)e* a; ^^ and a(k, t) = a(/c)e~ tw ^ fc ^. 
Putting it back into the expansion of $ gives rise to 

/oo 
dfccos(/b;)(a + (£;)e Mfc) * + a(k) e - iuj{k)t ) . (3) 
-oo 

The Fock Hilbert space 7i can be built up by acting by the creation operators on the vacuum: 

a(k)\0)=0 ; VA; 

\h,k2,...,k n ) = a + (h)a + (k 2 )...a + (k n )\0) ; h > k 2 > . . . > k n > . 

Note that in labelling the states, k is always positive.^] For technical reasons in some formulas 
we also allow k to take negative values, but we always mean a symmetric extension, that is 
a{k) = a(—k). The vacuum expectation value of the time ordered product 

/°° d 2 k j P ~ ik o(t~t') 
7^T2 >7^( ettl(M } + elkl(X+X } ) ( 4 ) 

solves the inhomogeneous equation: 

(□ + m 2 ) (0\T($(x,t)$(x',t'))\0) = -i8 N (x,x')5(t-t') . 

Besides the usual bulk propagator, which describes how the field propagates from (x, t) to (x', t'), 
(H) also contains another contribution, which can be interpreted as a bulk propagation of the 
field from (— x, t) to (x',t'). Thus the free boundary theory (Neumann boundary condition) can 
be realized by the mirror trick: We compute every quantity in the usual bulk theory, but any 
time we insert a field at (x,t) we insert the same type of field also at the mirror point (— x, t). 
Since the interacting theory is defined in terms of the free quantities, (in the calculations we use 
the free propagator) we have the following interpretational consequence: The particles interact 
not only with themselves but also with their mirror partners. 

1 actually k provides only a different parametrisation of the energy by the relation k = \Joj 2 — m 2 , since in the 
presence of a boundary the momentum is not conserved. 
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Interacting theory, asymptotic states 

Non trivial interaction is described by (0) when £/($), V($) or both are non zero. To handle this 
case we use the adiabatic hypothesis. That is the interaction is switched on adiabatically in the 
remote past and switched off in the remote future. Moreover, we also suppose that the particle 
spectrum does not change during this adiabatic procedure: Only the masses are renormalised. 
In a real scattering experiment the prepared state is one of the free theory (in state) and the 
detected state is also a free state (out state). Both the in and the out states provide a basis for 
the Hilbert space H, and the .R-matrix is a unitary transformation connecting the two: 

\final) out = R\initial) in . 

The unitarity of the i?-matrix expresses the fact that the transition probabilities sum up to one. 

The simplest process we can imagine is in which a single particle travels toward the boundary, 
reflects on it and then returns as a multi-particle configuration. The naive bulk analogue of this 
process is trivial since asymptotic particles are stable by assumption. In the presence of the 
boundary we can interpret this process in terms of the mirror transformation. In this language 
the in state contains not only the incoming particle but also its mirror image w.r.t. the boundary 
x = 0. The scattering of the incoming particle on the boundary has a contribution describing 
its scattering on its mirror image, a process analogous to the two particle scattering in the bulk. 
To be concrete: In the initial state of this process we have a wave-packet of the form 

/oo 
dk'f(k')\k'}i n , (5) 
-oo 

where we suppose that f(k') is well localised around k. The space time dependence of such a 
configuration is 

/oo 
dk'f(k') cos(k'x)e- iuj{k ' )t . 
■oo 

It describes a wave packet travelling with momentum k towards the boundary. It also contains 
however, the mirror image of the packet which is on the other side of the wall (so is not in the 
real space-time) and travels with momentum —k as shown on the figure: 




If there is no interaction (free case) then this state (f|) is the eigenstate of the free Hamiltonian. 
Since the time evolution is trivial the picture in the remote future looks like 
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Now the real and reflected particle travels away from the boundary with momentum | — k\ 
and the mirror image with momentum k. 

In the interacting case the final state may contain particles (or just one particle in the inte- 
grable case) travelling backward from the boundary. This coincides with the idea of (l| where 
the in state contains a particle with rapidity 9 while the out state with rapidity —9. 

Boundary reduction formula 

We are interested in the case when both the in and out state contain a single particle with 
definite energy. The energy conservation can be factored out: 

in {k'\R\k) in = 2n(5(k - k') + 6(k + k'))uj(k)n(\k\) . 

Our aim is to make correspondence with the other definitions of the reflection matrix. For this 
reason we express the reflection matrix 7Z in terms of the correlation functions. In the bulk 



theory this is done using the reduction formula ||14|| . In the following we derive an analogous 
formula for boundary theories. The steps of the derivation are similar to the ones in jTJ] . 

The in field is a free field so we can use the decomposition ([3|). The inverse of this relation is 

r° « 
a in (k) = 2i dxcos(kx)e luj(k)t dt^in(x,t) (6) 

J -co 

r° <-> 
a+{k) = -2i dxcos(kx)e- iuj(k)t dt®in(x,t) . 

J — oo 

Using the definition of the in state we have 

out(Pl,---,Pk\qi,---,qi)in =■ () =out (Ph ■ ■ ■ ,Pk\atn(ll)\<l2, ■ ■ ■ , qi)in ■ (7) 

Now apply formulas (H) to obtain 

f° -■ )^ 

() = -2i / dxcosiq^e luJ ^ 1 > t g t out (p 1 , . . . ,p k \$ in (x, t)\q 2 , ■ ■ ■ , qi) in ■ 

J — oo 

We suppose that the in field can be expressed in terms of the interacting field as t) — > 

Z 1 / 2 $j n (x,t) as t — >• — oo. As a consequence 

{) = -i lim Z- 1/2 2 / dxcosiq^e-^'dtoutiPi, ■ ■ ■ ,PkMx,t)\q 2 , . . . , gi ) m . 

t — > — CO I 

J —CO 

Since lim 3>(x,t) = lim Z 1 l 2 § out {x 1 i) we also have 

t^co t^co 

() = out (pi,---,Pfc|a+ t (gi)|g 2 , •••,%)« + 

/0 poo ^ 

dx / dtdticosiq^e-^^dtautipi, ■ ■ ■ ,PkMx,t)\q 2 , ...,©)<„} , 
-OO J —CO 

from which the connected part is 

^•0 rco 

iZ- 1,2 2 j dx I dte~ Mqi)t {cos( qi x)(out\d 2 t <&{x,t)\in) + {out\^{x,t)\in){-d 2 x + m 2 ) cos(gix)} , 
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where (out\, (\in)) is the shorthand form for (pi, . . . ,pk\, (\q2, ■ ■ ■ ,9/)), respectively. Performing 
the partial integration, (which is legitimate if momenta are smeared with wave packets of the 
form of d||)), we have to be careful to keep the surface term. The connected part turns out to be 

%Z- 1/2 2 ! d 2 xe- iuj( - qi)t cos( qi x){D + m 2 + 6{x)d x }(out\$(x,t)\in) , 



where / d 2 x = J°. dx dt is the integral over the entire physical space-time. This is the 
first stage of the reduction formula. In the second step we eliminate an outgoing particle. The 
derivation straightforwardly follows the combination of the previous computation and the usual 
bulk derivation. The connected part of the result is 



out(pi,...,Pk\qi,...,qi) t n = - 4 ^ _1 / d 2 xd 2 x , e i ^ t '-^^co S (q 1 x)cos(p 1 x) 

{□ + m 2 + 5(x)d x } {□' + m 2 + 6(x')d x > } (p 2 , • • • ,p„|T($(x, t)$(x , t))\in) . 

Iterating the steps above the general matrix element ([?]) can be expressed in terms of the k + I 
point function. 

In particular for the reflection matrix we have 

out < k'\k > in - out < k'\k >out= 2n(5(k - k') + 6(k + k'))uj{k){K{\k\) - 1) = (8) 
-4Z-/^Ve-M W ,'^,,,{ D + m , + S(x)d J s{a + m » + «(,^} G (,,,', t - , 

where 

G(x, x',t-t') = (0|T($(x, t)$(x, t'))\0) . (9) 
Perturbation theory 

Let us turn to the description of the interacting theory as a perturbation of the free one. In 
doing so we use the interaction representation. That is the time evolution operator is given by 
the time ordered (T) product as 

U(t) = Texp\-i [ dtH mt (t 

This Hamiltonian contains the in fields and acts on the in Hilbert space by construction. Clearly 
the i?-matrix can be expressed as 



R = [/(oo) = Texp |-«y dtH int (tyj 



which also gives a direct calculation of this quantity. The interacting field is built up from the 
free field as 

<$>{x 1 t) = U- l {t)§ m {x,t)U{t) . 
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Putting this expression into the two-point function (H) and using the usual heuristic derivation 
we obtain 

, , , , s ,< 'xx, v (0\T(<$> in (x,t)<& in (x',t')exp{i f d 2 xC int \<& in (x,t)] })|0) 

(0|T(exp{? J d 2 a;£ iri i[$ iri (x,t)J))|0) 

Making an expansion of the exponential we arrive at the perturbative series 

En=o£(°l r (Ms,t)$in(z , ,0 / rf 2 gi£ inf [^ n (gl,ti)].../d 2 X n £ iwt [^in(a:n,U])|0) , . 

In computing the vacuum expectation values of the product of the fields we can use Wick's 
theorem. The results are encoded in the Feynman rules which are given in Appendix A. 

From careful analysis of the perturbative series (see Appendix B) one can deduce that the 
momentum space Green function has the following form: 

G(p,p',uj) = 2n(S(p + p) + 5(p - p'))G bu i k {p,u) + G bu i k (p,u)B(p,p' ,u)G bu i k (p' \u) . (11) 

Here G bu i k (p,u) is the propagator of the bulk theory, which in terms of the spectral function 
a(m 2 ) has the usual Kallen-Lehmann representation ||T4]| : 

^ / x iZ f 00 , /2 ia(m' 2 ) . . 

G bu i k {p,uj) = — r-— + / dm — — — . (12) 

uj z — p z — m z + it J 4m 2 oj a — p z — m' z + le 

We also have the decomposition 

B(p,p',uj) = B 1 (p,p',uj) + B 2 (p,u) + B 2 (p',u) + B 3 (u). (13) 



The interpretation of the terms in ( |TT|) is the following: The first term describes the propagation 
in the presence of the boundary without hitting the boundary. In the second term G bu ik(p,u) is 
the propagator to the boundary, B(p,p',u) is the reflection on the boundary, while G bu i k (p',uj) 
describes the propagation back from the boundary. In the reflection matrix Bx(p,p',uj) really 
depends on both momenta and comes from the purely bulk interactions, B 3 (uj) is the purely 
boundary contribution and B 2 represents the cross terms. 

Now we are able to relate the two different definitions of the i?-matrix. Performing both 



inverse Fourier transformations of ( |TT|) in the momentum variables, but keeping only the contri- 



butions of the poles of the first term in propagators 0l2[) : 



G(x, x', t-t')= I - ^tV ^-*' 1 + (1 + ^( fc H;^)^) )e -^)(^)) , 
2n 2k(uj) 2k(u) 



where k(u) = \/uj 2 — m 2 . Comparing the reflected wave with the unreflected one the reflection 
matrix was defined to be 

= l + ZB{k{u) k(u),ui) 
2k{LU) 

We will recover the same result from our boundary reduction formula (||) . First we recall that 
the reduction formula describes the way how the matrix elements are related to the correlation 
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functions. Considering the correlation functions in momentum space the operator □ + m 2 gives 
a factor of — k 2 + m 2 for each external leg. The space-time integrations, as inverse Fourier 
transformations, put all the momenta on shell. Since all of the outer legs in the correlation 
functions are dressed up in the perturbation theory to contain the exact bulk propagators ( |T^ ) 
with poles of the form pz^r, the reduction formula merely amputates the legs and gives the 
residue of this multi-pole pole. In the boundary case we have an analogous interpretation. 
Similarly to the bulk case the momentum conserving part of does not give any contribution 
to the i?-matrix so it is enough to consider the other term. A careful analysis shows that 
(□ + m 2 ) is the operator which amputates the legs starting with a bulk vertex, while S(x)d x is 
responsible for amputation of the legs starting with a boundary vertex. As a consequence the 
(□ + m 2 )(D' + m 2 ) term gives Bi(k, k',u(k)), the terms (□ + m 2 )5(x')d x > and (□' + m 2 )S(x)d x 
together give B2(k,u(k)) + B 2 (k' ,u(k)), finally S(x)d x 5(x')d x > gives B 3 (u(k)). We also have an 
overall factor 2iiZS(u(k) —u(k')) expressing energy conservation. Collecting all these terms and 
using the identity 2n8{uj{k) — uj(k')) = ^p-2n(5(k — k') + 5(k + k')), we obtain that 



K(k) = 1 



k 

ZB(k,k,u(k)) 



2k 



which shows that the reflection factor defined with use of the asymptotic states and the one 
defined using the two-point function are identical. 



Conclusion 

The boundary reduction formula, derived in the paper for a boson with the most general (possibly 
non-integrable) boundary and bulk self-interaction, showed the equivalence of the previously used 
definitions for the i?-matrix, and supports the correctness of the formalism. This formulation 
enables one to derive the main properties of the i?-matrix such as analyticity, unitarity, crossing 
symmetry and analyse its pole structure directly without referring to the crossed channel picture 
used in The analysis of the perturbative series, Landau equations, Cutkosky rules, the 
derivation of the boundary Coleman-Thun mechanism and of the analyticity properties of the 
i?-matrix are the subjects of our next paper. 
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Appendix A 

The Feynman rules of the boundary theory can be obtained by expanding the interaction La- 
grangean as 

OO OO ft 

t) = E *)" + ^ E w m t)M (14) 

N=0 ' M=0 
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and putting flip into (|l0|) . In contrast to the perturbative approaches developed so far we 
formulate the resulting Feynman rules in momentum space. That is we compute 

/oo roo /»oo /*oo 

dx x . . . I dx n dh... dt n e ij: J^ x ^ tj) G(x h h . . . x n , t n ) 
-oo J —oo J —oo J —oo 

according to the Feynman rules: 

• Draw all possible oriented graphs with n external legs such that the external legs point 
inside the graph. The lines (both external and internal) can be of two types: straight and 
dashed. The vertices can be two types: black and white. Label the external legs by the 
2-momenta (pi, uji) . . . (p n , u> n ), and the internal lines by new momentum variables (k , k°). 



To a black vertex with iV incident lines attache the contribution 

str dsh 



ia N 2n5(Y, *? - £ k°)n5(J2 k\ k]) (15) 



i£in idout i£in i£in idout 



To a white vertex with M incident lines attach the contribution 

^ M 27r5($> - J>°) (16) 



i&n i£out 



P - 




M 



To a line of any type labelled by (k 1 , k°) attache the contribution 

i 



k 2 — m 2 + ie 



If the line is internal then integrate over k: J (f^- 
Sum over all topologically distinct diagrams. 



(17) 
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The notation i G in/out means that the line labelled by (kj,k^) is in-coming/out-going at 
the given vertex. The label str/dsh on the sum means that the summation goes only over 
straight/dashed lines. 

These Feynman rules need some explanation. The derivation is fairly standard: One derives 
the graph rules first in coordinate space. This is done by writing G(x\,ti . . .x n ,t n ) as a pertur- 
bative series. Each term of the series is represented by a graph by associating the various parts 
of the contribution to the vertices and lines of the graph. In particular, there is a space-time 
integration associated to each vertex. According to the two terms in (PI), there are two kinds 
of vertices: black (bulk) ones corresponding to the first term and white (boundary) ones corre- 
sponding to the second term. Owing to the presence of the factor S(x) in the second term of ( PI ) 
at the boundary vertices the space integration can be performed. There is only a single type of 
graph lines at this stage. 

Converting to momentum space (Fourier transformation of G(xx,t\ . . . x n , t n )), all the remain- 
ing space-time integrations can be performed resulting S functions at each vertex on the momenta 
of the lines incident to the vertex. At the boundary vertex there is only a time integration left, 
hence we obtain a 5 function involving only the time-like component of the momenta. At the 
bulk vertex there are both time and space integrations resulting in two 5 functions. The presence 
of both Xi — Xj and Xi + Xj dependent terms in the free propagator ([|) implies that we have a 
sum of 5 functions on the space-like component of the momenta. The 5's in this sum differ in 
the signs of some momenta in the argument. The two types of lines are introduced in order to 
associate different graphs to the individual terms: a line is straight if it corresponds to the term 
coming from the X{ — Xj dependent part of the free propagator associated to the given line, and 
dashed if it comes from the Xi + Xj dependent part. 



Appendix B 

In this appendix we analyse the properties of the two point function in momentum space. A 
systematic investigation of the correlation functions can be achieved by generalising the para- 
metric representation ((li]] section 6.2.3 ) to the boundary case. Since the introduction of all the 
machinery is quite lengthy and we need the result only for the two point function, instead of the 
detailed presentation we summarise how the formula ([Tl| ) can be obtained by the direct analysis 
of the perturbative series. 

Energy is conserved at each vertex so we have 

G(p 1 ,u 1 ,p 2 ,uj2) = 2ix5(uji - u 2 )G(p 1 ,p2,uj 1 ) . 

Momentum is not conserved in general. There are graphs whose contribution spoils momentum 
conservation: such is any graph containing a boundary vertex. On the other hand a graph 
with only bulk vertices and straight lines gives rise to a contribution proportional to S(pi — p2)- 
Furthermore, the contribution is proportional to the contribution of the same graph computed 
in the bulk theory. This can be verified by comparing the appropriate graph rules with the bulk 
theory rules, and one finds that the only difference is a factor 1/2 in the bulk vertex contribution. 
That is the contribution of a graph with Nf, bulk vertices, no boundary vertices and only straight 
lines is 2~ Nb times the contribution of the same graph was in the bulk theory. 

On the other hand there is a symmetry operation on the graphs that leaves the contribution 
invariant. Namely, at each bulk vertex we can change the types of all incident lines and also 
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the signs of the momenta labelling the out-going lines. One checks that this transformation - 
performed independently at each of the Nb bulk vertices - does not affect the graph contribution 
and gives rise to a symmetry factor 2 Nb . This symmetry factor just compensates the factor 2~ Nb . 
Hence the sum of the contributions of the graphs with only bulk vertices and straight lines plus 
the ones related to them by the symmetry transformation described gives the propagator of the 
bulk theory. It is not difficult to see that no other graph gives contribution respecting momentum 
conservation. 

Summing then the contributions of the graphs with only bulk vertices and straight internal 
lines (and arbitrary external ones) plus the graphs related to them by symmetry transformation 
we obtain the "momentum-preserving" part of the propagator 

2tt (5{pi + p 2 ) +6(pt - p 2 ))G bu i k (p 1 ,uj 1 ). 

Now we concentrate on the momentum non-preserving part. A diagram which remains con- 
nected when any of its internal line is cut is called one particle irreducible. Any diagram with two 
external legs can be built by attaching one particle irreducible diagrams after each other. For the 
momentum non-preserving diagrams we can separate the consecutive momentum preserving one 
particle irreducible subdiagrams adjacent to a given external line. This gives subgraphs that are 
identical to the graphs in the series of the bulk propagator, so they give factors Gb u ik{pi,^i) and 
Gbuik(P2, ^2) at the external lines, respectively. The contributions of the remaining momentum 
non-preserving subgraphs are collected in B(pi,p 2 , u>i). As a consequence we have the following 
form for the propagator: 

G(pi,P3,wi) = 2n(5(p 1 +p 2 ) + 5(pi -p2))G bu i k (p 1 ,Ui) + G bu i k (p 1 ,w 1 )B(p 1 ,p 2 ,u 1 )G bu i k (p 2 ,uj 1 ) . 

If in a Feynman graph, contributing to B(pi, p 2} uJi), both external lines are incident to a bound- 
ary vertex, then it does not depend on any of the momenta. Its contribution is collected in 
B 3 (u). Terms depending on one momentum only are collected in B 2 (pi, u)\), these are the graphs 
in which one of the external lines ends in a boundary vertex. The contribution of a diagram 
starting with bulk vertices on both ends depends both on p\ and p 2 . Such contributions are 
collected in Bi(pi,p 2 , cji). Summing up all these terms we have flT3|). 
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